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ABSTRACT 

Motivated by recent observations of solar hard X-ray bursts with a 
spatial resolution of a few arc seconds , we compute the theoretically pre- 
dicted spatial variation of hard X-ray flux versus height in both thick 
target non-thermal, and thermal, models of solar hard X-ray bursts. Our 
work on the thick target model revises previous results in this area by 
adopting a more realistic model for the flaring atmospheric structure, and 
also by taking into account energy loss and scattering processes in the 
evolution of the non-thermal electron beam which have previously been 
neglected. 

It is pointed out that in the so-called "thermal" model currently in 
vogue, there is a substantial non-thermal bremsstrahlung component at all 
photon energies; it is further shown that this non-thermal component results 
in a hard X-ray flux versus height distribution whose maximum brightness 
per unit length is, similar to the thick target model, located in the 
chromosphere, and not in the thermal source itself (where the bulk of the 
total low photon energy bremsstrahlung is emitted). The characteristics cf 
the hard X-ray height distributions in both thick target and thermal models 
are summarized and compared, and predictions as to the observable spatial 
structure of the X-ray flare in both cases are made. 


Subject Headings : Sun: chromosphere - Sun: corona - Sun: flares - Sun: X-rays 
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I. INTRODUCTION 

Current theoretical models of the hard X-ray bursts commonly associated 
with solar flares fall into two main classes - thermal and non- thermal. In 
the former model the X-rays are attributed to collisional bremsstrahlung 
within an ensemble of thermally relaxed electrons (e.g. Chubb, Kreplin, and 
Friedman 1966; Brown 1974; Crannell at *1 . 1978; Matzler g_t «1. 1978; Brown, 
Melrose and Spicer 1979; Smith and Lilliequist 1979; Emslie and Brown 1980; 
Emslie and Vlahos 1980; Smith and Auer 1980; Emslie 1981a). In the latter 
the X-rays similarly result fiom bremsstrahlung, although this time from a 
beam of suprathermal electrons incident on a cold background target. The 
exact target geometry involved has led to a variety of non-thermal models — 
thick target (Brown 1971; Emslie 1980), thin target (Datlowe and Lin 1973), 
and trap (Takakura and Kai 1966; Brown and Hoyng 1975; Melrose and Brown 
1976; Hudson 1978; Emslie, McCaig, and Brown 1979). 

Current X-ray spectral measurements (e.g. Hoyng, Brown and van Beek 
1976) are unable to distinguish between the various models outlined above, 
because of the instability of the X-ray spectrum deconvolution problem (e.g. 
Craig 1979) and because of the ability of a thermal distribution with a non- 
uniform temperature to mimic a "non-thermal" (e.g., power-law) spectral form 
(Brown 1974). Similarly, available polarization measurements (Tindo et al. 
j 970, 1972a, b; Tindo, Mandel'stam, and Shuryghin 1973; Tindo, Shuryghin, 
and Steffen 1976) are not reliable enough (Mandel'stam 1980) to discriminate 
amongst the predictions of various models (Brown 1972; Langer and Petrosian 
1977; Bai and Ramaty 1978; Emslie and Brown 1980; Leach and Petrosian 1981). 
It has therefore been suggested (e.g. Rust and Emslie 1979; Emslie and Rust 
1980) that spatially resolved hard X-ray images may provide vital clues to 
aid in the discrimination amongst candidate models. 
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Observationally, the only "spatially resolved" hard X-ray data published 
at the time of writing is the fortuitous observation of a behind-the-liob 
flare by two satellites separated by a suitable distance in ecliptic longi- 
tude (Kane e£ al . 1979). However, with the successful launch of the NASA 
Solar Maximum Mission satellite, with its Hard X-Ray Imaging Spectrometer (van 
Beek et_ al . 1980) on board, there is now an abundance of spatially resolved 
(£8" resolution) data available for study, and preliminary results have been 
presented (de Jager et al. 1981; Hoyng e£ al . 1981). 

Theoretically, the ~viy published study of the height distribution of 
solar hard X-rays in the context of a particular model is by Brown and McCly- 
mont (1975), who analyzed the thick and thin target non-thermal models. How- 
ever, the predicted distribution in these models is a sensitive function of 
the density structure of the model atmosphere adopted; Brown and McClymont 
(1975), in making use of a quiet Sun model (Vernazza, Avrett, and Loeser 
1973), in fact obtain results which significantly differ from those obtained 
by using a more realistic flare model atmosphere (Machado and Linsky 1975; 
Machado et al. 1980); we show this below in §111. 

In the present paper we therefore revise Brown and McClymont' s (1975) 
calculations, using a more realistic model atmosphere as a target, and also 
taking into account the effects of collisional scattering and reverse current 
ohmic energy losses on the non-thermal electron beam (effects neglected by 
Brown and McClymont in their analytic calculations) . We also compute the 
height distribution of hard X-rays in the dissipative thermal model (Brown, 
Melrose, and Spicer 1979; Smith and Lilliequist 1979; Smith and Auer 1980), 
giving due recognition to the non-thermal high energy tail component in the 
electron velocity distribution In such models (Brown, Craig, and Karpen 1980; 
Emslie and Vlahos 1980) . The hard X-ray height distribution in coronal trap 
models (Takakura and Kai 1966; Melrose and Brown 1976) depends on the magnetic 
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field geometry in the trap and on the electron scattering process within the 
trap (Melrose and Brown 1976; Wentzel 1976), and will not be dealt with here. 

In ill we describe the method by which our (numerical) results are obtained. 

In §111 we compare our thick target results with those of Brown and McClymont 
(1975). In §IV we study the thermal model and show that, contrary to some 
expectations (Rust and Emslle 1979; Emslie and Rust 1980), the hard X-rays 
in this model do not have their brightest component in the thermal part of 
the flaring loop, although the bulk of the total photon flux at low energies 
emanates from the (larger) thermal region .In §V we discuss the results obtained 
and show how the two models may be distinguished by spatial measurements with 
sufficiently high resolution; such spatial resolution is within current instru- 
mental capabilities. 
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II. METHOD OF COMPUTATION 

Since we shall assume the hard X-ray emission from any thermal source 
to be spatially uniform (see §IV), it remains only to describe the computa- 
tion of the non-thermal component present In both the thick target and 
dissipative thermal models. We shall make use of the procedure developed 
by Emslle (1978, 1980) to compute the dynamics of, and X-ray bremsstrahlung 
yield from, the beam of non-thermal _*ectrons; we note that his mean scattering 
and energy loss theory yields results which are very similar to a full Fokker- 
Planck analysis of the beam-target interaction (Leach and Petrosian 1981). 

It has been pointed out by Chambe and Henoux (1979) that Emslle' s 
(1978) treatment of ion-neutral collisions neglects the inelastic momentum 
loss from the charged beam particles. This loss is given by 



Inelastic 


1 dE 
mv tl dt ’ 


( 1 ) 


where V|| is the component of velocity in the precollision direction, and 
m and E the beam particle mass and energy respectively. Adding these In- 
elastic momentum loss terms to Emslle 's (1978) equation (22) gives the 
revised total momentum loss rate in a target of arbitrary Ionization level: 



dt 




(1-x) 



( 2 ) 


In this equation m^ and m^ are the electron and proton masses respective- 
ly, e (e.s.u.) is the electronic charge, x is the fractional ionization 
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number of hydrogen, and A, A', A" are Coulomb logarithms, defined in Fmslie's 
(1978) paper. It follows that Emslie's (1978) equation (25) for the scatter- 
ing parameters 8 should read 


« - 2xA + (1-x) A” 

S e A* + x (A-A*) 



(3) 


(8 and 8 corresponding to electron and proton bombardment respectively), 
e p 

In turn, equation (14) of Emslie (1980), which includes the effect of both 
collisions and reverse current ohmic energy and momentum losses from the 
beam electrons, should read 


&Lm 

dN 


-te't 


2xA + (1-x) A 


"] 


+ 



(4) 


where \i is the direction cosine (w.r.t. the vertical) of the beam electrons, 

n and n are the plasma density (cm”*) and resistivity (e.s.u.) respect- 
* 

ively, E q is the minimum injection energy required of an electron to pene- 
trate to a column depth N (cm 2 ), and F , E and <5 define the injected 
electron flux spectrum: 

F (E ) (electrons cnf^s^keV -1 ) - (6-1) F E 6 " 1 (E +E )“ 5 . (5) 

o o oo oo oo o 

Note that Emslie's (1980) equation (13) remains unaltered. 1 
*We have re-performed Emslie's (1980) calculations using the revised version 
of his equations (13) and (14) and find no significant changes from his results. 

Using the revised equations as outlined above, we may now derive the 
variation of electron flux with depth in the target, and so compute the X-ray 
bremsstrahlung yield as a function of depth and photon energy (see Emslie's 
[1980] equation [22]). Emslie (1S80) integrated this over depth to obtain 
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Che resulting bremsstrahlung spectrum only; however, It is a straightforward 
matter to use his results directly to compute the height variation of the 
X-ray flux at a given photon energy. Results follow in the next two sections. 
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III. RESULTS FOR THE THICK TARGET MODEL 
As mentioned In $1, th« earlier treatment of this problem by Brown and 
McClymont (1975) uses the quiet Sun atmospheric density structure (Vemazza, 
Avrett, and Loeser 1973) as a target. This structure In fact bears little 
resemblance to that found In flares (Machado and Linsky 1975; Litas and Cook 
1979; Machado et al . 1980). The principal difference between the quiet Sun 
and flare model atmospheres Is the particle column depth (cm z ) of the 
transition zone, with its associated steep density gradient; In empirical 
flare model atmospheres N T2 a 10 J Vr 2 , while N^ is *2 x 10 l8 cm 2 In the quiet 
Sun (Vemazza, Avrett, and Loeser 1973) and =3xl0 19 cm 2 in active regions (Basri 
et al. 1979). It is not immediately obvious whether one should employ an 
active region or flare empirical model as a target for the bremsstrahlung- 
producing non-thermal electrons, since they are flash-phase phenomena and 
empirical flare modeling to date (e.g. Machado and Emslie 1979; Machado et al. 
1980) applies only to the later stages of the flare. However, since In a 
thick target model the principal cause of the increase in N^ z is the bombard- 
ment by the bremsstrahlung-produclng electrons, and since it Is easily shown 
that preflare chromospheric material rises in temperature by the required 
10 7 K or so in a very short time (see eq. [8] of Emslie and Noyes 1978), 
it seems more reasonable to use the flare models of Machado £t al . (1980) 
as a background target for our computations. (For further discussion of 
the validity of empirical flare model atmospheres, see Emslie, Brown, and 
Machado 1981.) Now, the predicted height distribution of thick target hard 
X-rays is a combination of two effects: (i) density variation (the thick 

target yield is proportional to the ambient number density), and (ii) 
electron flux attenuation (the yield is also proportional to the electron 
flux at a given level, which is in turn a decreasing function of the column 
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density). Thus ons sees that by moving tht lsvsl at which the larga photon 
flux anhanca&ant occurs dua to affact (i) to a largar (flara) valua of the 
column density N (and so to a laval at which tha electron flux has attenuated 
to a much higher degree), a significant difference from Brown and McClyaont * s 
(1975) predicted X-ray flux versus column density profile may result. 

In addition to these considerations, we note that Brown and McClymont 
(1975) neglect colllslonal scattering of the beam electrons, and also 
reverse current ohmic enorgy losses from the beam; these processes can 
result in substantial modifications to the electron flux versus depth profile 
(Emslie 1980). Further, Brown and McClymont* s (1975) assumption of a fully 
ionized target underestimates the bremsstrahlung yield at greater depths, 
where the fractional ionization level of hydrogen is small (Brown 1973). 

All these effects are incorporated in the present treatment. 

We have calculated the variations of hard X-ray flux versus depth for 
two model atmospheres. The first uses the model chromosphere 'FI* of 
Machado et el. (1980) and the second the model chromosphere 'F2' of these 
authors (hereafter these models will be referred to as MAVN1 and MAVN2 
respectively). To these models we attached isothermal coronae (temperature 
3 x 10 7 K) 3 x 10*km in length (measured down one half of the loop). The density 
at the loop apex was taken to be 3 x 10 lo cm * for model MAVN1 and 10 ll cm”* 
for model MAVN2, and the density variation along the coronal part of the 
loop was taken to be smoothly Increasing with depth such that the integrated 
column density at the transition zone matched Machado et al_. *s (1980) 

values of 1.5xl0 20 cm 2 and 1.6 x 10 21 cm 2 respectively. These coronal 
models are admittedly somewhat crude; however, it turns out that the coronal 
hard X-ray emission per unit length in the thick target model is negligible compared 
to that in the chromosphere (Figures 1 and 2), so that any deficiencies in the 
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the model adopted will not significantly affect the conclusions of the paper. 

In Figures 1 end 2 va show the variation of the fraction of the total 

X-ray flux emitted per km of vertical distance at depth z (km) (c ilumn depth 

H fc.-’] ), for these two flare models. The electrons were assuawd to be 

Injected by a "black box" mechanism (with no intrinsic hard X-ray signature) 

at the top of the coronal loop (N4) with a flux spectrum given by equation 

(5) with F • 10 1 * electrons caT*e~* # E *20 keV, and S ■ 4 (see Emslle 
oo oo 

1980). The dashed line indicates a purely colllslonal calculation of the 
electron beam dynamics and the solid line a calculation involving reverse current 
losses as well. These reverse current losses are quite sensitive to the atmospheric 
model adopted, particularly in the corona (Emslle 1980), which is vhy we show 
the purely colllslonal calculation for comparison, even though It Is 
unphysical. Note also that the requirement that the reverse current be 
stable to current-driven plasma turbule sets an upper limit to the 
significance of reverse current ohmic dissipation in the electron beam 
energetics. Only for certain coronal conditions (e.g. high density, approx- 
imately equal electron and ion temperatures ) is auch a large flux F qq 
allowed to pass stably; if these conditions are not met then F qo must be 
reduced and the effects of reverse currents correspondingly decreased 
(Emslle 1981b). Note finally that the figures only extend down to a depth N = 

10 22 cm” 2 . This is the thermalization depth of electrons with initial (i.e. 
Injected) energy %300 keV (e.g. Emslle 1978); since data on hard X-rays 
with photon energies of this order are poor, we consider extrapolation to 
higher electron energies (and so greater depths) than this highly uncertain. 

If the spectrum (5) does in face extend to arbitrarily large energies Z q , 
then <t>(N) will asymptotically approach a power-lew form with index (1-6) 12 
• -1.5 (Brown and KcClymont 1975; Emails 1978). 
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Figure 2 


It la evident from the figures that the effect of reverse current 
ohmic dissipation Is significant for the coronal region of model MAVN1 
only; this Is due to the relatively low coronal densities and consequently 
small colllslonal losses (compared to ohmic ores) In this region (Emails 
IV 80) . Mote that the effect of the reverse current Is to enhance the 
coronal component of the hard X-ray emission relative to the chromospheric 
component. In fact (see Emslle 1980). both chromospheric and coronal 
components are reduced (In absolute flux units) from their purely colllslonal 
values; however the effect is more severe on the chromospheric emission due 
to the reduction In flux penetrating to large depths effected by the reverse 
current ohmic dissipation in the corona. 

Comparing our results to those of Brown and McClymont (1975) we notice 
that in model MAVN2 (Figure 2) the region of brightest X-ray emission per 
unit length Is narrower than In either Brown and McClymont' s (1975) results 
or the results for model MAVN1 (Figure 1) and is located at the base of 
the transition region, as opposed to mid-way Into the chromosphere. This 
is because, as remarked above, the transition zone level in mcdel MAVN2 
is so deep that the electron beam has already significantly attenuated 
before it reaches the chromosphere, to such an extent that the increase of 
chromospheric density with depth from the base of the transition zone down- 
wards cannot overcome this diminishing flux effect. In the quiet Sun and 
MAVN1 models, however, this coronal attenuation is not so significant, and 
consequently the region of maximum bremsstrahlung yield per unit length is 
pushed downwards into the chromosphere and made broader. 

These differences are, however, probably quite undistinguishable by 
observation, at the current level of instrument capability. The height 
range over which $ ^ 0.1$ is =2xl0 5 km (3" in angle) for Brown and 

DIAX 
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McClymont's quiet Sun calculation (see their Figure 4), =*10 3 km (IV) for 

model HAVN1 (Figure 1) , and *5xl0 2 km (1”) for model MAVN2 (Figure 2) . All 
of these angular sizes are below the currently available spatial resolution 
of hard X-ray detectors (van Beek et al . 1980) . 

It is worth commenting further on the results for model MAVN2. Figure 2 
shows that the depth dependence of hard X-ray emission in the coronal part 
of the loop is significantly different for low (e =»• 10 keV) and high (e s 70 
keV) photon energies. For low energies the principal cause of the 
decrease of <P with depth is effect (ii) above, namely electron flux attenuation, 
while for high energies effect (i) (increasing density with depth) dominates, 
causing the X-ray flux per unit length to increase with depth. One would expect 
effect (ii) to begin to dominate effect (i) (causing a maximum in <{0 when one 
reaches the depth at which electrons of energy E K e begin to lose energy rapidly 
(note the dependence of the Bethe-fleitler bremsstrahlung cross-section [Heitler 
1954^ on electron energy E) . This occurs when (e [keV] ) 2 £ 6* e 4 AN * 

(N/10 17 ) (Emslie 1978); the curve for e = 40 keV in Figure 2 has a maximum 
at N = 3 x lO 20 cm 2 , in good agreement with this prediction. 

This concludes our discussion of the thick target model. In the next 
section we shall consider the thermal model, and compare the results from 


the two models in § V. 
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IV. RESULTS FOR THE THERMAL MODEL 

In this section wa investigate the height distribution of hard X-rays 
expected in the "dissipative thermal" model, originally proposed by Brown, 
Melrose, and Spicer (1979), and currently under extensive investigation (see 
Smith and Lilliequist 1979; Ensile and Brown 1980; Emslle and Vlahos 1980; 

Smith and Auer 1980; Smith and Brown 1980; Emslle 1981a). This model is 
characterized by a "thermal" distribution of electrons at the top of a coronal 
loop, confined parallel to the loop axis by a strong magnetic field and per- 
pendicular to the loop axis by a pair of collisionless ion-acoustic turbulent 
fronts, which result from the instability of the reverse current associated 
with the outward heat flax from the source (see, e.g. Kahler 1975) and 
which move along the arch at approximately the local ion-sound speed (Smith 
and Lilliequist 1979; Smith and Auer 1980). 

It was recognized by Brown, Melrose, and Spicer (1979) and by Vlahos 

and Papadopuulos (1979) that in such a model very high energy electrons are 

not confined by the turbulent plasma fronts and escape downwards in much the 

same manner as a thick target beam. Although the exact details of this escape 

process are somewhat uncertain (Smith and Brown 1980), it is typically found 

that all electrons with velocities v £ 3 v (where v is the mean electron 

e e 

thermal velocity) can escape. Brown, Craig, and Karpen (1980) and Emslie and 
Vlahos (1980) have considered the thick target bremsstrahlung yield from 
these high energy electrons and find that it is comparable to the thermal yield 
for photon energies e greater than a few times (the mean electron 

thermal energy). Specifically, Emslie. and Vlahos (1980) give the ratio of 
thick target to thermal bremsstrahlung at photon energy t as 


n(e) 


i 6.6 X 10' 1 8“ 


nL 4(6-2) 


e ~8 2 /2 e e/kT (kT) 2 ln 




; e < E 


min 


( 6 ) 


; e > E 


min 
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igure 3 


igure 4 


where k is Boltzmann's constant, B the ratio of the critical escape 

velocity to v and 6 the spectral Index of the power-law escaping tail 

(see Emslle and Vlahos 1980); note that assuming the tall to be a simple 

extended Maxwellian [see Brown, Craig, and Karpen 1980] only reduces n by 

a factor of ^ 2 for the parameter set considered and so does not affect 

the conclusions to follow . n (cm -3 ), T(K), and L(km), are the density, electron 

temperature, and length of the thermal source and is the threshold energy 

for an electron to escape through the conduction fronts ■ *s8 z kT. We shall 

hereafter consider a source model with 8*3, 6*4,n* 10 11 cm 3 (see Smith 

and Lilllequist 1979; § III) and T • 22.2 keV (so that the electron escape 

energy threshold E . is 100 keV) ; results for other values of these parameters 
min 

can readily be obtained using equations (6) and (8). For the parameter set 
chosen. 


n(e) 


369 

L 


exp (0.045 e)ln(^~) 


(7) 


(with e measured in keV) ; this is shown graphically in Figure 3. 

If we assume that the thermal bremsstrahlung emission is spatially uniform 

over z < z* (N < N*) , then the fractional (relative to the combined thermal plus 
thick target non-thermal emission in one leg of the arch) yield per km of loop is 

! 2 + n(e)4> (e,z)L 

LL . z < z* 

[l + n(e)]L 

( 8 ) 

n Cc) 4> (e,z) 

Li ; z > z* 

1 + n(e) 

where $ (£,z) is the thick target fractional yield as found in Sill, but 

with a lower cutoff of 100 keV introduced to the spectrum (5). 

In Figure 4 we show the resulting <Xzl structures for various e; results 
for the thick target model neglecting reverse currents (see Figure 2) are shown 
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for comparison. Model MAVN2 was used In Che calculations; the results using 
model MAVN1 are similar. The column density N* was taken to be 10 2 ° cm 
and the corresponding source length 2*10 cm (note the symmetry in the 
source model) . Note that the ratio of thick target to thermal values of 
$ in the model is in fact Independent of L : to see this we rewrite equa- 
tion ( 7 ) as 

n 0 Ce) 

n(e) • — r~ » (9) 


with 

n 0 (e) - 369 exp (0.045 e)ln(^~) * 


( 10 ) 


Then equation ( 8 ) becomes 


<Me,z) 


EC 


L + n (e) 

O 


; z < 


n 0 (e)* tt (e,z) 
L + n Q (e) 


Z > Z* 


(ID 


so that the ratio of <p values in the 'throoospheric"and "thermal” parts of 
the source is 


9(e) 


n o (em tt ^ £> ^chromosphere 
2 + thermal r.gion 


( 12 ) 


Since for all values of e the second term in the denominator is negligible 
compared to unity (the suprathermal electrons give off little bremsstrahlung 
in the relatively tenuous coronal plasma - note the sharp drop off in at 
z> z* - Figure 4), we see that 0(e) is, to a good approximation, constant : 



16 


0(e) 


V e)[ *tt ^’^chromosphere 
2 


(13) 


It Is immediately evident from Figure 4 that the maximum bremsstrahlung 
yield per unit length at all energies still occurs at the base of the transi- 
tion region . This is contrary to earlier speculations (e.g. Rust and Emslie 
1979; Emslie and Rust 1980) that the most intense hard X-ray emission should 
be emitted in the thermal region at the top of the arch. Comparison of the 
thick target and thermal model results in Figf e 4 shows that the fractional 
contribution from the coronal region is indeed enhanced in the thermal model; 
however this is not at the expense of the chromospheric emission but rather 

at the expense of the emission in the region between the conduction fronts 

2 

and the transition region (i.e. between . » -2 x 1(T km and z ■ 0) . 

"Note that <j> (z) is an increasing function of z for all energies in this 
region of the atmosphere; only when N reaches =10 21 cm 2 , the thermalization 
depth for electrons of energy E 4 - 100 keV, does any significant decrease in <p 

occur; cf. discussion of thick target model curves in § III. 

This behavior is in fact easily understood: since the electron energy required 

to penetrate to chromospheric depths in model MAVN2 is ^ 80 keV (see, e.g., 
Emslie 1978), the fractional chromospheric yield in the thick target model is 
matched in the thermal model by the contribution from the high energy escaping 
tail. Further, due to the high background plasma density in the chromosphere, 
the high energy tail thermalizes in a very short distance and also has a high 
bremsstrahlung yield per unit length (the additional fact that the chromo- 
spheric plasma is only partially ionized further enhances the chromospheric 
emission - see Brown 1973). However, in the tenuous coronal plasma between 
the confining turbulent conduction fronts and the transition region, the 
thick target component, being inefficient compared to thermal emission at 
these densities (Smith and Lilliequist 1979; Smith and Auer 1980) is very 
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small. We thus s«e that the height distribution of hard X-rays in the dissi- 
pative thermal model ia characterized by two very low intensity regions, one 
in each leg of the flaring loop, flanked above and below by relatively bright 
regions, making five regions in all. The more Intense of these bright regions 
are the lower ones, with the degree of contrast between them and the uppermost 
region increasing with increasing photon energy (see Figure 4). 
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V. DISCUSSION 

The principal result of the above analysis is that the brightest hard 
X-ray region is always at the base of the transition zone, in both thick 
target and thermal models. This result follows in the thick target model 
because of the use of a more realistic model atmosphere than previously 
employed (cf. the results of Brown and McClymont 1975), and in the thermal 
model because of the presence of a precipitating high energy non-thermal tail 
of electrons, which behave exactly as a thick target beam. 

In our discussion of the thermal model (§IV) we have tacitly assumed 
that the high energy tail is always present (i.e. we have assumed a steady- 
state situation). Processes by which this escaping tail may be repopulated 
from the bulk thermal plasma are not well understood (seeEmslie and Vlahos 1980; 
Smith and Brown 1980) and so it must be conceded that perhaps the X-ray height 
distribution depicted in Figure 4 strictly applies only to an initial tran- 
sient state. In such a case the "relaxed" <f> distribution would simply be a 
rectangle of length L/2 and height 2/L (see Figure 4). Clearly whether or not the 
high energy tail is continuously repopulated is of crucial importance; for 
the present we note that the interpretation of microwave bursts by Emslie 
and Vlahos (1980) in the context of the dissipative thermal model demands 
that the tail does exist. 

Returning to the results of Figure 4, we note that the distribution will 
in fact be time-dependent, as the conduction fronts advance along the arch 
(recall, however, the results of § IV [ e^uativ : (12) ], 
which showed that this variation of L does not affect the dominance of 
the chromospheric component in the $[z] profile). An obvious corollary 
to this is that the source has only a finite lifetime T ^ L^/c^ , where 
L^ z is the distance from the energy release point to the transition zone 
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L -1 

and c la the ion sound speed * 10* T* cm a . For plausible parameters 
s 

we see that T Is of order 30 seconds, which Is somewhat greater than the 
time resolution of available hard X-ray imaging Instruments (van Beelc et al. 
1980). Thus it appears that the introduction of a "multiple kernel model", 
such as has been advocated by Brown, Craig, and Karpen (1980) to explain 
observed (rfitzler at al. 1978) emission measure/temperature correlations (see, 
however, Emslle 1981a) may be unnecessary. (Note, however, that even if such a 
multiple kernel model were appropriate, with the kernels all formed near the 
top of the arch, we should replace L by a mean L in the analysis of §IV, 
thereby resulting in a distribution very similar to that of Figure 4 [recall 
the insensitivity of the results to L — § XV ] ) . 

Finally we note that the best spatial resolution currently available 
in the hard X-ray spectral region is rather limited (^, 8" ; van Beek et al. 

1980) . An angular resolution of this order corresponds to a length of 

= 5000 km on the solar surface, implying that we should really scan Figures 

1, 2, and 4 with a rather broad "filter" of this width. The effect of this 
smoothing depends on whether we are considering thermal or non-thermal 
models, and on the model atmosphere adopted. For instance, for the thick 
target model with atmosphere MAVN1 (Figure 1), one coronal "pixel" will 
contain some 10~ 5 x 5 x 10 3 = 5% of the emission, while a chromospheric 
pixel will contain some 70X of the total emission (from one leg of the arch). 
Thus in this case we should expect to see two chromospheric "bright points" 
at the feet of the flare loop, with little emission elsewhere. 

In the thick target model with atmosphere MAVN2, we find that, due to the 
greater amount of coronal material compared to atmosphere MAVN1, a coronal 
pixel now observes some 15X of the emission (at £ - 40 keV), while the 

chromospheric pixel observes a similar fraction of the total emission. Thus 

we should here expect to see almost umtorm emission over the whole loop. 
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Finally, In the thermal model, we see that a pixel observing the thermal 
region will contain some 30-40* of the emission, a chromospheric pixel will 
observe some 20-30*, and a pixel in the region between the thermal region and 
the transition zone will observe at most 2 or 3 percent of the emission. In 
this case we therefore expect to see the characteristic "five layer" structure 
mentioned in I IV. 

These observational predictions are sufficiently distinct to distinguish 
between the above three scenarios with currently available instrumentation 
(van Beek et al. 1980), and it is honed that results from the Solar Maximum 
Mission Hard X-Ray Imaging Spectrometer will be able to cast some light upon 
which scenario actually applies to solar flares. 


I thank J. Leach and P.A. Sturrock for their comments on the manuscript 
and help in clarifying the presentation, and the referee for his excellent 
criticism of the original manuscript. This research was supported by NASA 
Grant NGL 05-020-272 and ONR Contract N00014-75-C-0673 . 
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FIGURE CAPTIONS 

Figure 1: Fraction of total thick targat bremsstrahlung yield $ (km* 1 ) 

varaua depth t (kai) and coluan depth N (cm” 2 ) for photon anargiea 
10, 40, and 70 keV with modal atmosphere MAVN1 (Machado at al . 1980) 
as background. The z-coordlnate la measured from the transition 
zona downwards; note tha dlffarant scala in the coronal and 
chromospheric regions of tha figure. Tha solid curves include 
tha affects of both collislonal and reverse currant ohmic loasas 
from tha breasstrahlung-producing non- thermal electron beam; 
the dashed line takes into account only collislonal losses. 

Compare with Figure 2 and with Figure 3 of Brown and McClymont 
(1975). 

Figure 2: As for Figure 1, except using atmospheric model MAVN2 (Machado 

et al. 1980) as background. Since MAVN2 has a deeper transition 
zone than MAVN1 the electron flux has already been significantly 
attenuated by the time it reaches the transition zone; thus the 
region of maximum bremsstrahlung yield per unit length is very 
localized and occurs at the base of the transition zone. This 
result is markedly different to the results of Brown and McClymont 
(1975), due to their adoption of an unrealistic (quiet Sun) 
background atmosphere. 

Figure 3: Ratio of thick target to thermal hard X-ray bremsstrahlung (n) 

as a function of photon energy c in the "thermal" model dis- 
cussed in §1.'. Above the electron escape energy threshold of 
100 keV there is Insignificant thermal emission and n becomes 
infinite. Note the substantial contribution of the thick target 
emission at all energies; this has important implications for 
the height structure of hard X-rays in this model (see §IV and 
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Figure A). 

Figure 4: <t> versus N for thermal model with parameter set outlined in 5 IV. 

The figure is normalized so that jw* is unity for one leg of 
the arch. The hard X-ray flux versus height structure is 
characterized by two dark regions (one in each leg of the arch), 
flanked outside and in between by bright regions. Note that the 
4> (N) structure Is time dependent due to the motion of the conduc- 
tion fronts through the corona, although this does not affect the 
relative sizes of 0 values in the corona and chromosphere (see 
§ IV) . 
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